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A tensor network formalism of thermofield dynamics is introduced. The formalism relates the 
original Hilbert space with its tilde space by a product of two copies of a tensor network. Then, 
their interface becomes an event horizon, and the logarithm of the tensor rank corresponds to the 
black hole entropy. Eventually, multiscale entanglement renormalization anzats (MERA) reproduces 
an AdS black hole at finite temperature. Our finding shows rich functionalities of MERA as efficient 
graphical representation of AdS/CFT correspondence. 
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I. INTRODUCTION 

Applications of anti-de Sitter space / conformal field 
theory (AdS/CFT) correspondence [l[ to statistical and 
condensed matter physics are hot topics in string theory. 
On the other hand, it is recognized that a new class of 
variational anzats in statistical physics, so called mul- 
tiscale entanglement renormalization anzats (MERA), 
would be a discrete version of the AdS / CFT correspon- 
dence 0-0| ■ Both of them overcomes difficulties of real- 
space renormalization in critical systems. It is thus at- 
tracting attention to examine their complementarity in a 
mathematical level. 

The MERA is a kind of tensor network (product) 
states (TNS, TPS) of quantum many-body systems on 
lattices |7H10| . Historically, the TNS formalism was con- 
structed so that the ground-state variational wave func- 
tion satisfies the entanglement-entropy scaling. The ex- 
amination was first developed for gapped cases. In the 
gapped cases, the entropy obeys the well-known area- 
law scaling [Ih-[l5j. In spatially one dimension (ID), 
the wave function appropriate for the scaling is matrix 
product state (MPS). Actually, the MPS is numerically 
optimized by the density matrix renormalization group 
(DMRG) method, and the DMRG is known to be the 
most powerful method in quantum ID systems [IB - ll8j ]. 
A natural generalizartion of the MPS to higher dimen- 
sions is to make the tensor contraction that represents 
a set of short-range entangled pairs on bonds. In that 
sense, the TNS is also called projected entangled pair 
state (PEPS). When the tensor rank and the surround- 
ing area of a partial system are m and A respectively, the 
entanglement entropy between the partial system and its 
enviromnent is given by See = Alnm oc A. On the 
other hand, in critical cases, the area-law scaling violates 
logarithmically fl9l - l29j |. Then, we must take a large m 



value of the order of the total sites L, if we keep the ba- 
sic TNS structure. However, it is possible to construct a 
hierarchical tensor network in one higher dimension that 
has a managable size of the tensor rank and shows the 
logarithmic divergence of the entropy. That is a concept 
of the MERA network. 
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FIG. 1: 2D binary MERA network. Filled dots, diamonds, 
and filled triangles are the original sites in a quantum ID 
critical system, disentangler tensors, and isometries, respec- 
tivery. The vertical direction denotes renormalization flow. A 
red line represents the causal cone enclosing a partial system. 
The surface area of the causal cone is given by the sum of the 
number of the boundary points in each layer. 

The MERA network has its graphical representation as 
shown in Fig. [I] and it is easy to visualize the entangle- 
ment structure inherent in the network. A key ingredient 
of MERA is that the network is composed of layers, each 
of which has its own energy and length scales. The lay- 
ered structure emerges as a result of repeating the block- 
spin transformation as well as disentangling transforma- 
tion. In an another viewpoint, the tensors with large \ in 
the TNS is decomposed into a set of tensors with smaller 
dimensions and with different functionalities. Then, the 
total entanglement entropy is given by the sum of the en- 
tropy in each layer compatible with the area-law scaling. 
The reason for the appearance of the area law is that the 
smaller tensor rank represents more classical-like short- 
range correlation in this extended space. Graphically, the 
sum corresponds to the outside area of the causal cone in 
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the discrete AdS space, and actually this is comparable 
to the Ryu-Takayanagi formula 30] . 



II. TENSOR NETWORK REPRESENTATION 
OF THERMAL STATE 



In this paper, we focus on general TNS and in particu- 
lar MERA at finite temperature. In comparison with the 
ground states, finite temperature properties of the TNS 
are less understood. Because of the potential complemen- 
tarity of MERA and AdS/CFT, the understanding also 
gives us deeper insight for application of the AdS/CFT 
to condensed matter physics [3ll - l34j . In the application, 
it is necessary to deform the asymptotically AdS metric 
so that the IR geometry has a black hole solution. The 
black hole is a source of coarse graining, and determines 
the temperature at the conformal boundary of the AdS 
space. Then, we would like to know whether the black 
hole naturally emerges from the TNS formalizm without 
any phenomenological assumptions. This is the purpose 
of this paper. 

In a viewpoint of the CFT, the TNS has rich function- 
alities. Let us briefly look at a basic property of the MPS 
that is the most primitive TNS in ID. The property is 
about a dicrete version of the Calabrese-Cardy formula 
for the entanglement entropy [2(| [2l[ , the so called scal- 
ing of entanglement support given by See = (c/c/6) lnx, 
where c is the central charge, K — Q/c{y/l2/c + 1) is the 
finite-entanglement scaling exponent, and x is the matrix 
dimension [35l - l37| . It has been shown that the correla- 
tion length £ is given by £ = \ K . This means that the 
MPS is a very simple structure but catches the essential 
feature of the CFT. Therefore, we expect that the black 
hole is also described by the TNS. We will confirm this 
statement with use of the MERA network. 

Since the MERA is a wave-function approach, it is 
straightforward to use thermofield dynamics (TFD) for 
finite-temperature formalism of quantum field theory 
(the readers may be aware of equivalence between TFD 
and dual CFT in a special case) [38h43| . For our purpose, 
we reformulate the original TFD so that the TFD wave 
function matches well with the tensor network represen- 
tation. We will find that by this reformulation the event 
horizon naturally appears as a result of hidden quantum 
entanglement between the original Hilbert space and its 
tilde space. The entanglement produces a combined ten- 
sor network, and we will examine the basic properties of 
the network. 

The paper is organized as follows. Section II is the 
main part of this paper. We first examine a single-site 
model, and find a method for constructing the TNS for- 
malism of the thermal state. Then, we examine the 
finite-temperature MERA network. We will present mi- 
croscopic derivation of the black hole entropy, the event 
horizon, and the temperature at the conformal boundary 
of the AdS space. In Sec. Ill and Sec. IV, we discuss 
related topics and summarize our study. 



A. Entanglement between Dual Hilbert Spaces and 
Thermal Entropy 

Let us start with a particular model in TFD. In addi- 
tion to the original Hilbert space, we introduce the tilde 
space that is isomorphic and decoupled to the original 
space. Then, it is able to represent a thermal average by 
the expectation value of the thermal state. We define a 
generator G and its evolution operator U(6) by 



G 
U{6) 



iftbj (aa — a' a') , 

i8G 



(1) 
(2) 



where the parameters = (3/2 = l/2fcsT and ui charac- 
terize inverse temperature and the energy scale of exci- 
tation modes, respectively. The Fermion operators a and 
a satisfy the anti-commutation relations 



Here, we assume 



{a, at} = {a, a/ \ = 1, 
{a, a} = {a, a} = 0. 

{a, a} = {a, a' } =0. 



(3) 
(4) 

(5) 



In the following, we define the vacuum state by |00) = 
|0) <g) |6>, at |0) = 0, and at |6) = 0. We use the notation 

(6) 
(7) 
(8) 



a f 00) = 


|io), 


at 00) = 


01> , 


tat 00) = 


|11). 



Let us consider the following thermal state 
\O(0)) = U{6) |00) 

= (u(9) + v{9)c$tf) 

where we can find 



1 00), 



u(6) 
v{6) 



cos (Hui8) . 
sin (fiwff) . 



(9) 
(10) 



(11) 
(12) 



In order to examine the entanglement entropy of this 
system, we define the partial density matrix p by 



p = fr\0(9))(0(6)\ 
= u(6) 2 \Q) (Q\ + v(0) 2 



1)1 



(13) 
(14) 



where Tr traces over degrees of freedom in the tilde space. 
In a context of TFD, this is nothing but the thermal den- 
sity matrix. At the same time, p represents the amount of 
entanglement between the original and the tilde Hilbert 
spaces. The entanglement entropy is given by 



See = -Tr(plnp) 

= -u{8) 2 \nu{ef 
< In 2. 



v{e) 2 \nv{ef 



(15) 
(16) 
(17) 
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On the other hand, the thermal entropy of this system is 
clearly (the system takes |0) or |1)) 



St = kB In 2. 



(18) 



Identifying S EE with S T /k B [14|, \M IMtti, we find 
hujQ = 7r/4 and 



T 



2huj 
Trk R ' 



(19) 



which maximizes See an d then See looks like St- 



entangled. Futhermore, there is a fact that the entropy of 
extremal black holes known so far has been explained by 
assuming maximally entangled states, and this fact also 
matches with our situation. In that sense, our entropy 
seems to behave as the black hole entropy and the event 
horizon is located at the interface of A m and A n . Then, 
Eq. (fl9|) corresponds to the Hawking temperature of our 
system. The maximally entangled system appears in the 
equal probability case (pi — 1/x for i = 1,2, ...,%) even 
when the degree of freedom increases. Then, the upper 
bound of the entanglement entropy is given by 



B. Vector Product Form of Thermal State 

In order to examine the physical meaning of Eq. (fTJ 
we first transform Eq. (jlOp into a vector product form. 
That is represented as 



|0(6>)> = J2 AmAh \ mfl ) 



(20) 



m,n— 0, 1 



where A m and A n are x-dimensional vector and its tilde 
conjugate, respectively. The index m (h) takes or 1. 
Usually it is considered that the original and the tilde 
Hilbert spaces are independent. However, this hidden 
correlation mediated by the rank \ plays an important 
role on the emergence of the black hole. The tilde conju- 
gate of the vector is not defined in the original TFD, but 
in this paper we define it by the transposition as well as 
the complex conjugate of the vector elements. This wave 
function is nothing but the MPS. Here we consider the 
single-site problem with open boundary condition, and 
then the edge matrix is terminated by an appropriate 
vector in order to get a scalar coefficient. 
We introduce the following representation 



( A i\ 



A m = (A?,A!?,--- ,A™),A n 



A', 



(21) 



Then, we obtain the following solutions 

A°=( v ^),o),A 1 = (o, v / ^) 



and 



A 



o 







(22) 



(23) 



We can find more general solutions, but according to 
symmetry of two Hilbert spaces, we take the vector ele- 
ments so that A n becomes a copy of A m . In the present 
case, the proper (minimal) \ value is 2. When we com- 
pare Eq. ([14]) with Eqs. ([22]) and ([23]), we clearly find that 
the thermal entropy counts the degrees of freedom at the 
interface between A m and A n when they are maximally 



EE 



i=l 



Pilnpi = lnx, 



(24) 



and this also leads to the thermal entropy. 
In general, the thermal state is represented as 



= p 1/2 \i) 



(25) 



where p is the density matrix, and we take \I) 
according to the general representation theorem [4£ 
decompose p into singular values {A,,} as 

p 1/2 = VAV, 



(26) 



where V and V are column unitary matrices, A is a diag- 
onal matrix and diagA = (Ai, A2, ...). Then the \ value 
is given by 



X 



= rank(UA 1/2 ) = rank(A 1/2 y). 



(27) 



C. 



Decomposition of Vector into Truncated MERA 
Network and Emergence of AdS Black Hole 

For general L-sites systems, Eq. ([20]) is extended as 

ini2---rn L j^ht n 2 — fit 



i*> = £ E Amir 

{rrij } {fij } 

x \mim 2 ■ ■ ■ m L hih 2 ■ ■ ■ n L ) 
where rrij and fij are local valuables and 



(28) 



a=l 



It is argued in ref. [50j, [5l| that truncating some tensors 
from the complete MERA network roughly represents 
the AdS black hole. We derive it from a microscopic 
viewpoint. In our notation, A™ 1 ™ 12 '" mL corresponds to 
the truncated tensor network. Thus, we decompose it 
into a set of tensors with smaller dimensions. The de- 
composition of A™ im2 '" mL into any network is always 
possible. Depending on criticality of our target models, 
we select an appropriate network. The tensor geometry 
should match with the symmetry of the original quantum 
system. 
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FIG. 2: Decomposition of the vector A into a set of tensors. 



Hereafter we focus on the conformal field theory de- 
scribed by the 2D MERA network. Then, we decom- 
pose ^™i™2 - mt mt0 a MERA-like network. It is noted 
that the presence of the index a leads to the incomplete 
MERA network. Figure [5] is an example that is the four- 
sites MERA network under the open boundary condition 
but the top tensor is truncated (the network is termi- 
nated). This diagram is mathematically represented as 



.4 



m 2 ,m 3 



W 1 - TJ™ 2m3 W 5 

r r m\m2 m2m3 r y m^m^i 



(30) 



with use of isometry W and disentangler U. Since the 
vector rank is %, the indices 7 and 5 take 1, 2, y/x- 

We generalize the above result. The total degrees of 
freedom at the interface is originally but let us assume 
that the tensor rank and the number of the tensors at the 
interface are respectively m and A after the decomposi- 
tion. Then, x is represented as 



X = m 

Here, A is related with L by the following condition 
L 



(31) 



A, 



(32) 



where th is the layer number starting from 0. In the 
binary MERA network, we take 77 = 2. In general, each 
renormalization process merges 77-sites together. Substi- 
tuting Eqs. (|3"Tj) and (|3"2")l into the upper bound of the 
entanglement entropy, we obtain 



'EE 



lnm. 



(33) 



Now our discrete AdS space in Fig. [T] is represented by 
the metric 



ds 2 



{d(r In 77)} + (77 T dx) 2 , 



(34) 



and changing the valuable r to z = rf (z is normalized 
by the AdS curvature), we obtain the standard notation 



9 dz 2 + dx 2 
ds 2 = - . 

z z 

Thus, Eq. Q33p is transformed into 



'EE 



lnm, 



ZH 



(35) 



(36) 



where zh = rf H . According to the CFT at finite tem- 
perature, the entanglement entropy of the original ID 
quantum systems is given by 



See = 77 In [ — sinh ( ^ 
3 \TTe \ p 



(37) 



where e is a UV cutoff. Assuming that (3 is small enough, 
we expand Eq. (|3T[) as 



(38) 



The second term is proportional to the system size L, 
and thus See obeys the volume law at high temperature. 
When we identify Eq. (J36} with Eq. ([551. we find that 



(39) 



fes T=(-ln™Ucx: 

.C7T I Z H 



-1 
H ■ 



This is nothing but the temperature scale arizing from 
the AdS black hole, and is consistent with the pre vious 
results based on the holo grap hic principle (3ll434l l52j . 
Our case is similar to ref. [53| in a sense that the black 
hole entropy is interpreted as the entanglement entropy 
in CFT living on the boundary of the AdS space. In our 
formulation, the separation of the CFT and the event 
horizon was derived from decomposition of A™ im2 "' mL 
into the truncated MERA network. 

Before going to the next subsection, we further com- 
ment on the physical meaning of Eq. (I39|) . We require 
consistency between our result and the holographic the- 



ories up to the coefficient, 
should be given by 



Then, the coefficient of 



3 1 1 
— In m = — , 

C7T 27T 



and this requires that 



m 



c/6 



1. 



(40) 



(41) 



Here we have supposed the minimal series (0 < c < 1) 
and the Gaussian CFT (c = 1). The right hand side of 
Eq. (|40[) has been obtained for continuous systems. Thus, 
Eq. (|41|) is just a criterion of the m value. However, the 
result indicates that the MERA generates a really 'clas- 
sical' space near the event horizon, since the quantum 
entanglement has been almost vanished for the m value. 
This small m value might be the reason for the success of 
semiclassical treatment of the black hole thermodynam- 
ics. 



D. MERA - Tilde MERA Combined Network 

Up to now, we have considered a method for decompo- 
sition of A™ im2 '" mL into a terminated MERA network. 
Let us imagine that the coefficient of the basis in Eq. (|2"8")l 
is a product of two terminated MERA networks. In the 
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FIG. 3: Combination of MERA (lower half) and tilde MERA 
(upper half) networks. A red wavy line represents an event 
horizon. Two pairs of isometries across the horizon are en- 
tangled. 




FIG. 4: Higher temperature case of MERA - tilde MERA 
combined network. When we fix system size L, larger \ leads 
to the smaller layer number. 



tilde-space side, we call it as 'tilde MERA' network. Fig- 
ure El shows schematic representation of our extended 
MERA network. We also show higher temeperature case 
in Fig. SI The total system is composed of the original 
and tilde MERA networks. At finite temperature, some 
upper layers of the original MERA network is truncated 
as already discussed. Two truncated MERA networks are 
pasted together at the position of the event horizon. The 
surface area of the horizon is determined by Eq. (1321) , and 
the temperature of the conformal boundary is given by 
Eq. (|39[) . In ref. [ilj . a similar viewgraph is introduced in 
a context of gravitational collapse in the AdS spacetime. 
Thus, our result is a strong support of similarity between 
MERA and AdS/CFT correspondece. 

The most recent development of the MERA is to con- 
struct the branching MERA network that reproduces 



anomalous entanglement-entropy scaling [54|. The net- 
work looks like asymptotically the AdS space, but the 
IR geometry is separated into various types of multiple 
branches. In this case, the global structure of the MERA 
- tilde MERA combined network is topologically different 
from the present case. It will be an intersting work to 
examine how the genus affects the excitation modes on 
the conformal boundary. 



III. DISCUSSION 

Based on the present results, we give some insights for 
black hole thermodynamics and for numerical renormal- 
ization group in condensed matter physics studied with 
use of TFD or dual CFT. 

It has been pointed out that the Kruskal transforma- 
tion in general relativity is similar to the Bogoliubov one 
in TFD [39H42T|. Both of them use hyperbolic functions 
for the basis transformation. When we look at the Pen- 
rose diagram, the diagram separates the whole spacetime 
into four blocks. They are two subspaces, and black and 
white holes. In between two subspaces, the black hole 
appears. This seems to be quite similar to our combined 
network. When an observer is in one of two subspaces, 
it is considered that the observer can not contact with 
the other subspace. However, we think that the observer 
indirectly contacts with the other space through hidden 
quantum entanglement. This information corresponds to 
the black hole entropy. 

The TFD has been also applied to the finite- 
temperature DMRG in quantum spin chains 55]. The 
natural ordering of sites for the best DMRG performance 
is to set up site 1, ancilla 1, site 2, ancilla 2, etc. The 
DMRG is a kind of diagonalization of the Hamiltonian 
with use of the reduced basis set. The basis set is deter- 
mined by the eigenvectors with large eigenvalues of the 
reduced density matrix that is obtained by tracing over 
environmental degrees of freedom in the whole system. 
Then, we would like to reduce the quantum entangle- 
ment between near-neighbor sites as much as possible. 
Since the state n and its ancilla h are maximally en- 
tangled at high temperature, they should be paired in 
the simulation. This means that in the MPS optimized 
by the DMRG the ancilla blurs the pure-state proper- 
ties on each site. On the other hand, in our combined 
MERA network, the bulk structure of the MERA net- 
work is modified. In that sense, the role of the tilde 
space on the coarse graining by finite temperature looks 
quite different in classical description (MERA) and di- 
rect quantum description (MPS). This may be a kind of 
duality between them. 



IV. SUMMARY 

Summarizing, we have shown that the dual tensor net- 
work formalism by the TFD is efficient for describing the 
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finite-temperature properties of quantum systems. The 
TNS is a simple tensor product, and we do not assume 
any gravitational setup. However, the black hole thermo- 
dynamics is automatically emerged from the formalism. 
In particular, the MERA network would be really repre- 
senting the discrete AdS/CFT correspondence, and then 
it seems that Eq. (j52")l is a very strong constraint. In 
2D maximally symmetric spaces, the negative curvature 
space only matches with the CFT in a group theoretical 
viewpoint, since the space is scale-invariant. Even if our 
MERA space is discretized, we think that the matching 
still remains. In the previous theories associated with 
MERA and AdS/CFT, the renormalization flow mech- 



anism was mainly examined. On the other hand, the 
present study has approached the global spacetime struc- 
ture of our extended MERA network. Thus, we hope 
that their interaction will leads to better understanding 
of MERA and AdS/CFT. 
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